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The spontaneous breaking of chiral symmetry and the U(l) axial anomaly in QCD, can be phe- 
nomenologically understood by considering instantons as the gauge configurations mediating quark- 
quark interaction. The existence of an exact zero mode solution of the Dirac equation in the field of 
a single instanton is the fundamental ingredient. Explicit expressions for ipo are available for T ^ 
and fj, = 0, and /i 7^ and T — 0. In this paper we derive the solution for the most general case 
T 7^ and p 7^ 0. This opens the possibility of investigating the QCD dynamics associated with 
instantons in the full phase diagram. As a first step in this direction, we will study the dependence 
of the instanton density from the thermodynamic coordinates. 



I. INTRODUCTION 

Our present understanding of the strong interaction is limited by the lack of analytical tools to investigate QCD in 
the infrared regime, where the theory is non-perturbative. Presently, the only rigorous non-perturbative framework 
to solve QCD is represented by ab-initio lattice simulations; however, in this framework, it is difficult to gain sufficient 
insight about the nature of the gauge field fluctuations which dominate the dynamics. This is essentially due to the 
method generally adopted for generating the statistical ensemble, since Monte Carlo importance sampling algorithms 
does not yield the configuration details needed for calculating the expectation values of interest. 

To shed some light on the mechanisms that underlie the QCD non-perturbative regime, a plethora of different 
models have been developed, in which the dynamical effects obtained can be clearly connected with the presence of 
specific degrees of freedom in the vacuum. The instanton model is a paradigmatic example of such models. 

The instanton is a well known semiclassical gauge configuration with integer topological charge, corresponding, 
in Euclidean space, to an exact self-dual solution of the Yang-Mills equations of motion. In the seventies a deep 
connection was established between QCD instantons and important properties of the theory, like chiral symmetry 
breaking and the U{1) axial anomaly [Tl[2]. It turns out that both phenomena are strongly related to the existence of 
zero modes and near-zero modes of the Dirac operator. In particular the violation of the axial charge in an arbitrary 
gauge potential is given by Aq^ — 2Nf{nL — where rtL,fl is the number of left and right-handed zero modes of 
the Dirac operator. 

Since an (anti-)instanton has one and only one exact (right-)left-handed zero mode, the anomaly emerges naturally 
at the level of the single instanton calculation. On the other hand, chiral symmetry breaking can be related to a finite 
density of quasi-zero mode through the Bank-Casher formula (qq) = — 7rp(A = 0), where p is the density of modes 
with eigenvalue A. The distribution of near-zero mode predicted by the Instanton Liquid Model, where a dilute and 
globally topologically neutral ensemble of pseudoparticle is considered, perfectly agrees with the phenomenological 
expectation and with calculations in Chiral Perturbation Theory [IHS]- 

More in general, the model provides a successful framework in the phenomenological study of the chiral dynamics 
of QCD, as proved extensively in literature [51-fTT]. 

In this work we will investigate the QCD phase diagram in the instanton model. For achieving this, one needs to 
extend the latter model to allow the study of the thermodynamic quantities. This requires, not only the analytic 
form of the instanton field at finite temperature (T) [12] . but also the general expression for the zero mode at finite 
temperature and chemical potential (/i). To date, the zero mode is known analytically only for the case of finite 
temperature with vanishing quark chemical potential and finite quark chemical potential with zero temperature |13l - 
I19j : as a consequence the phenomenology of the instanton model has not been extended to the full phase diagram. One 
attempt in this direction was pursued by Shuryak and collaborators [18j who, by extrapolating across the full phase 
diagram the results obtained along the axis, tried to extract general information on the chiral symmetry restoration. 
Nevertheless the absence of an analytic expression for the zero mode in the case of both finite T and finite fj, prevented 
to clarify in an unambiguous way the prediction of the model. 

In this paper we present the analytic expression of the zero mode as obtained through solving the Dirac equation 
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in the background of a single instanton for the general case of T ^ and pq ^ 0. With this new ingredient, all the 
terms entering the partition function of the model are known, and the study of the phase diagram of the instanton 
ensemble can be in principle carried out. As a first step we compute here the dependency of the instanton density 
from the temperature and chemical potential. The relevance of this last analysis can be understood looking at the 
relation (qq) = — l/(7rp)(3rt/2)^/^ which connect the order parameter of chiral symmetry breaking (qq), with a finite 
instanton density in the vacuum n. 

The paper is organized as follows. Sec. |ll]is devoted to the analysis of instanton induced fermionic zero modes. 
A review of the main features of these modes is presented in Sec. |II A[ where we spend some time focusing on the 



property of fixed chirality of the zero modes. Next, in Sec. II B we present the main result of the paper: starting 
from an Ansatz for the fermionic field profile, we solve the Dirac equation in the single instanton background at finite 
temperature and chemical potential, showing that the proposed profile is indeed a zero mode. Next, in Sec. |III[ we 
use the result obtained in the first part to start proving a systematic study of the phase diagram of the (instanton) 
model. To begin with in Sec. Ill A we give an explanation of the different terms entering the partition function of the 
model, to show how the knowledge of the analytic expression of the zero mode makes it possible to study the system 
for any T and /i. In Sec. |IIIB|we then show the instanton density in the full phase diagram and connect it to chiral 



symmetry breaking. Finally we present our conclusions and outlooks in Sec. IV 



II. FERMIONIC ZERO MODES 



Instantons are self-dual classical solution of the Yang-Mills Euclidean equation of motions with fixed integer topo- 
logical charge. In this paper we concentrate on singular gauge instanton, though all the results can be simply mapped 
to the regular gauge case. 

The instanton field is written in the form 

Alix)^-fj;,d,\nU{x), (1) 

where ^JJ^ is the 't Hooft symbol (main properties of it being presented in Appendix |A]) while the choice of the profile 

function n(a:) is dictated by the requirement of self-duality of the field strength F^jj[x) = F^j^{x). Indeed, substituting 
(jlj) in the definition of F^^, one has 

F;A^) = nlpi-d^dpXnllix) + {d^\nTi{x)){dp\nIl{x))) 

-(/iO^/)+7?;,(apln^(x))^ (2) 

and self-duality is guaranteed only if it is such that 

□ In n{x) + {dp In n{x)f = 0. (3) 

The choice of 11(2;) = 1 + jx^ satisfies the condition (Is]) and the parameter p identifies the size of the instanton. 



A. Chirality of the zero mode 

This paper is mainly devoted to the study of the zero mode of the Dirac equation as induced by instantons, at 
finite temperature and chemical potential. In this section we start the analysis reviewing the proof of one of the 
main characteristic of this fermionic modes: zero modes, associated with the presence of single instantons in the QCD 
vacuum, have fixed chirality. This information will be helpful in the following, searching for the analytic form of the 
zero mode. 

An important and direct consequence of this property lies in the very specific form of the instanton mediate quark 
interaction: two quarks exchanging an instanton necessarily have to flip their chirality. This effect is a very peculiar 
feature of the instanton gauge field, and has allowed to prove the presence of a finite density of instanton in the ground 
state of the theory [5] . 

By definition a fermionic field '^^{x) is a zero mode if 

iJl)i,{x) = Z7^(a^ - iA^,{x))i,{x) = 0, (4) 

where A^j is the instanton profile previously introduced. To prove that such a zeromode would have fixed chirality, it is 
useful to start the analysis from the square of Dirac operator. Eq.|4]can be squared yielding the equation JfP"il){x) = 
where the square of the Dirac operator Jf)^ is given by 



{^p^D^ f = -r-^l^,l.\F^,{x). 



(5) 
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At this point the adoption of a quaternionic notation, presented in Appendix [B] results particularly convenient. In 
this context the commutator of the gamma matrices can be expressed in terms of the 't Hooft symbol and ct^ matrices 
as 



r 1 oY^M-^" 



Throughout the paper we discuss the case of massless quarks, where the Dirac spinor can be written as a sum of a 
right and left handed mode: ip — ipR + tpL- This means that for the definition of the squared Dirac equation we have 
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- —v'^^F^^" 




(6) 



Considering the field strength (j2| and using the properties of the 't Hooft symbol (Appendix [A|, the term a°'r]'j^^F'^'^ , 
which relates to the right handed component of the spinor, can be rewritten in a more compact form as 

= ^[□inn + (a,inn)2]. (7) 

Here we recognize the left hand side of the self-duality condition ([s]); since for an instanton □ InH + {dp InH)^ = 0, 
we have that for the right component of the quark field the equation 

DVi? = 0, (8) 

holds. Moreover, since D is an hermitian operator, its square is positive definite, so that the equation ([s]) has no 
solution. The conclusion is that, if the zero mode exists, it must be left handed and such that it satisfies the equation 



- -Y^l.F^.J V'L = (9) 
The fixed chirality of the zero mode will help us in the research of the analytic expression for the zero mode. 

B. Zero mode profile at T 7^ and /i 7^ 

In this section we present our main result, namely we derive the quark zero mode profile in the general case of 
non-vanishing temperature and chemical potential. The starting point is the expression of the zero mode for the two 
limiting case T = 0, 7^ and T 7^ 0, = g]: 

Mr,t) = 2^e-*v/llM)^ i^f'*) 
where 7± projects on the desired chirality, and 



1 



( -io-2 



is a Dirac and color spinor. The other two functions entering expression ( 10 ) are the n(x) instanton profile, which 
here differentiates in the case of zero and finite temperature: 

1 H T = 

r2 + 12 

n(r,t)-<( ^ , (11) 

Tip^ sinh(27rrr) 

1 + — ^ T > 

Tr cosh(27rrr) - cos(27rtT) 
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and the new profile function $(r, t) defined by 



cos(r/x) + - sin(r/i) yU 7^ 0, T = 



cos(7rtr) 
cosh(7rrT) 



(12) 



The rest of this section is devoted to prove that, in order to obtain the generahzed zero mode at non vanishing T 
and n, we only have to replace the <i>(r, t) above with the Ansatz 



/cos(7rtr) cos(r^) sin(7rir) sin(r^) \ 
, = (n(r, t) - 1) (^—^^^-^ + j 



(13) 



To achieve our goal we proceed by inspection, proving that ( [13^ is indeed a zero mode solution, i.e. that: 

~ ^lJ4Mx■,^l,T)^o. (14) 

It is important to notice here that, in the case of finite chemical potential, the solution of the adjoint Dirac equation, 

V-^ {x; ~fi, T)i{^- i^^ - M74) = 0. (15) 

carries the chemical potential argument with opposite sign. This is crucial in order to have a consistent definition of 
the expectation values at finite ^. In particular looking at the case of zero temperature this property renders finite 
the norm, 



d'^xip^ (x; —fi)'ip'' (x; fi) — 1, 



whereas without the extra sign one has 



d^x^/'^(x; /i)'0^(x; /i) = c». 



(16) 



(17) 



This singularity being connected with the BCS singularity one encounters when resumming an effective (attractive) 
particle-particle interaction around the Fermi surface. 

To begin with we pass to the quaternionic formalism and write the zero mode profile as: 



et'^y/n ( 



01 



13] 



ia2 



p 



'fit 



n 



.do 



Making then use of the properties of the Pauli matrices 

= -0-2 and crj = -cr2Crj^o-2, 
we obtain the explicit expression for the color (j) and spinor (/?) components of the zero mode, namely 



(18) 



(19) 



At this point, we can use the fact that the zero mode is left-handed; projecting on the left-hand component of (19) 
we find 



1-75 
2 



[0 4^2]^,3ap( -e 



n 
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We can now proceed to introduce this result into the Dirac equation, obtaining 



-ia 
iuo 



p 



al3 



+ ^ i^l^t^ - o-l<^t^)di, hi n ) - ^ 



-ia, 
ia4 



al3 



Using the relations ( B2||B3P we obtain 



(0 -i-o,,5p(|' 



-fit 



and 





ia4 



(0 ala,)^^ = { ala, )^.^ 



From this equalities we get the Dirac equation in the following form 



-IfJL -}^a\On 



a2 =0. 



A straightforward calculation shows that the first term of the equation (21) becomes 



+ 7(4^. - cjla^)d. Inn) f e^*Vn ^(Tta, 



" " I n 



9„n 



n 

while for the second term we have 

ip e^Vn ald^ 



-fit 



= i2^e^*Vna4($/ne-'^*) + iAie^*Vncr^(9^($/ne-''*). 
Neglecting the normalization factor {2^/2^Tp)^^ , the sum of the two terms yields: 



2a„n a„ ( 



-/it 



n □ f -e-^* 
n 



Performing the derivatives and using that Dn = 0, we finally obtain: 

n"^ (□$ - 2^94$ + ^2$) = 0. 



(20) 



(21) 



(22) 



(23) 



(24) 



(25) 



Substituting the expression of $(r, i) as given in (131 we immediately see that equation (251 is satisfied, which 
concludes our proof. 

It is interesting to notice that ( 25 ) corresponds to the massless Klein-Gordon equation at finite temperature and 
chemical potential. We also observe that ( 13 ) is nothing else but the Klein-Gordon propagator for anti-periodic 
boundary condition in the t direction, being such that <I>(i) — — $(i + 1/T). On the other hand, the Klein-Gordon 
propagator with periodic boundary conditions is given by 



$s(r,t) = (n(r,t)-l) cos(r^) + 



sin(7rtT) sin(r/i) \ 
sinh(7rrr) / ' 



(26) 
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and considering the limit for vanishing chemical potential we obtain nothing else than the 11 (r, t) used in the instanton 
calculus: 

^ , N , T^P^ sinh(27rrT) ^, . 

lim $5 r, t) = l + ^ — — ^ ^— — = n r, t) 27 

1 r cosh(27rrJ ) — cos(27rti ) 

We thus see that both the instanton profile and the associated fermionic zero mode, are written in terms of 
propagators of Klein-Gordon scalar particles with appropriate boundary conditions. 



III. PHASE DIAGRAM OF THE INTANTON LIQUID 



The knowledge of the analytic form of the zero mode opens the possibility to explore the full phase diagram of the 
instanton model. In this section we use the result obtained in the previous section and derive the partition function 
of the model showing that this can be numerically evaluated also at finite temperature and chemical potential. 

Next, we will present some results on the phase diagram of the model, and in particular the dependence of the 
instanton density from the temperature and quark chemical potential. This knowledge is important not only to 
understand the structure of the QCD vacuum, but also in connection to chiral symmetry breaking, since the instanton 
density appears in the mean field relation: 

(28) 

As a result, we will see that, as T and /i grow, instantons tend to organize in molecules and the density n appearing 
in the formula is not the total density of instanton, but only the fraction of instanton that does not coalesce. 

Let us stress from the outset that the study performed in this section is for illustrative purposes and no conclusion 
on the phase transitions in the full QCD phase diagram will be attempted. On the other hand these results can be 
seen as a first attempt to study the effects related to the topological structure of the QCD vacuum on the phase 
diagram. 



A. Instanton partition function and overlap matrix elements 



In the instanton model the quantum field theory describing the QCD vacuum is reduced to an effective theory 
where a finite number of effective degrees of freedom, the collective coordinates of the instantons fli = {zi^ pi,Ui} 
(position, size and color orientation), are treated in a statistical mechanical framework. In this context the partition 
function of the instanton ensemble is given by 

where the sum runs over the total number of instantons (A^+) and anti-instantons (-/V-). The single-instanton ampli- 
tude n{pi) contains the semiclassical tunneling rate (including one-loop quantum corrections), as well as the Jacobian 
arising from the introduction of collective coordinates |2j. The instanton interaction can be divided into a gluonic 
part Sint and a fermionic part represented by the determinant of the Dirac operator of the instanton liquid model. 

The fermionic determinant, in the chiral limit, can be calculated in the subspace of zero modes; by neglecting higher 
order corrections, in particular 

det{tp - */i74) « det ^^^^^^ ^^-^(q^' ^) ) , (30) 
where the Tia matrix elements are given by the overlap integrals 

nl/T 



[/A 



dr J dxTp\{T — Z4, X — z; fi,T){ip — ifx^ijipAiT, x; ii,T). (31) 

The nonhcrmiticity of the finite-/^ Dirac operator is reflected by the fact that 

TAi{fi,T) = Tj^i-fi) + tUp,T), (32) 
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i.e., the fermionic determinant is complex, entailing the well-known "sign" -problem in the partition function. 

The overlap matrix element depends on the distance between instantons, the instanton sizes pi^A and the relative 
color orientation U — ujUA, Which can be characterized using the four-vector Ui, = l/(2i)Tr([/ cr^). The definite 
chirality of the zero modes (in the limit of vanishing current quark masses) entails that the I — I and A — A matrix 
elements are zero. In vacuum, Lorentz invariance implies that the overlap matrix element can be characterized by a 
single scalar function (201 [5T], e.g. Tja = iu ■ zf{z). At T, fi ^ this is no longer true, and the matrix elements have 
the more general structure 

TiA = iu4fi +i{u- z)f2 (33) 

with fi = fi{\z\, Zi, pi,rhoA) and z — z/\z\. These functions can be calculated by writing the zero mode solutions, 
derived in the precedent section, in the form tjjf ~ (l)i,{x){'y,y)ijU'^^x'j- In the general case /i,T ^ 0, we have 

(l)^{x) = 5yi(t)i{r, Xi) + 5uifi(j)r{r, Xi). (34) 

and in terms of (j)nu, the overlap matrix elements are given by 

TiA = 2m„ j^'^ dr j Ax{cp^d,ct>i - ^^d^cj^i + cj^iO^ 



- So^iKK + C0i) (35) 

where = (f)^{x — z) and (pl = (puix). Comparing (33 1 with the above equation we find that the functions /i,2 can 
be expressed by the integrals: 

fi = 2j^^ drj dx(^^tidi(l)i + 2'l^i/r + dr(bi) 

f ^^{dr^i - d^cfi) - M - ^rr^ • f ^Ui) ] , (36) 



r - z\ 








+f ■ z(l)i{d4(t>i - dr<t>i) - M M • r<t>i<t>i + z ■ -(f>^(l,i . (37) 



For the subsequent study of the thermodynamics of the model we need to integrate the fi 2 functions in the Z4 and 
r coordinates, which means that for each value of temperature and chemical potential we need to have a sufficient 
number of points in the (24, r) plane to ensure integral convergence implying in turn large computational times. As 
an example we show in Fig. [l] the functions /i^2 evaluated at T = 65 Mev and /j. = MeV. 



B. Chiral phase transition at finite temperature and chemical potential 

In the instanton liquid model the study of QCD thermodynamics is reduced to a statistical mechanics treatment 
of an ensemble of classical pseudoparticles. In this context spontaneous chiral symmetry breaking in the vacuum 
is generated by randomly distributed uncorrelated instantons-anti-instantons which allow for a delocalization of the 
associated quark quasi-zero modes corresponding to the formation of a nonzero (qq) condensate state (see, e.g., @] 
for a review of the argument). This can be understood, by thinking in terms of quarks that can travel arbitrarily 
long distances by randomly jumping from one instanton to another (anti-instanton) and thus may carry their chiral 
charge to spatial infinity. 

In this picture chiral restoration can in principle proceed in two ways: either instantons disappear altogether or 
they rearrange into some finite clusters which no longer support any finite (qq) condensate. 

From the most recent calculations of the most refined version of the instanton model at finite temperature |22j . 
it emerges that instantons survives well beyond Tc, and the same looks plausible also for the case of finite chemical 
potential. This suggest that chiral restoration will appear as the consequence of the rearrangement of single instantons 
into I-A molecules p4l[T5l[23). 

The investigation of the interplay between the different components in the instanton partition function at finite 
temperature and chemical potential can be addressed, by adopting a simplified version of the instanton model, the so 
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/i [GeV] 




FIG. 1: /i (upper panel) and /2 (lower panel) functions evaluated as a function of 24 and r on a grid (blue points). The 
solid surface, obtained interpolating the points, is plotted to show that the grid is sufficiently dense to reproduce smoothly the 
details of the functions. 

called "cocktail model" [HIITS]. In this model, which we follow closely, the ensemble is decomposed into a mixture of 
random ("atomic") and "molecular" configurations, which yields a grand canonical partition function of the form 



E 



(38) 



In the thermodynamic limit V4 — > 00, and using the Stirling formula, the thermodynamic potential becomes 



\n[Z 



a+ml 



V4 



-n„ In 



eza 



rin In 



(39) 
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The atomic and molecular "activities" are 

= 2Cp'-^e-'"^^dp,rt{T)(^ J d^zdU[TjA{T,p)TAi{T,ti)]p^\ (40) 
zm = C^p"^'-^^e-^''^-^dp,rt{Tf{[TiA{T,ti)TAi{T,fi)ff) 

^ C^p'^'-^^e-^''-H,,rt{Tf j d''zdU[TIA{T,n)TAI{T,^J,)f'p'^^. (41) 

where dpert{T) is a perturbative suppression factor which has the form 

, , ^ f 2Nc + Nf, / Nf 

dpert{T) = exp ( '-^{npTf - h + - ^ 

(- ln(l + iTTpTf/3) + 0.15/(1 + 0.15(7rpT)-3/2)8)^ . (42) 

The underlying approximation in this approach is that the values of the hopping amplitude Tja in each individual 
configuration are replaced by a product of their mean square values in an uncorrelated ensemble. Another important 
approximation introduced in the model is the reduction of the gluonic interaction to an average repulsion term 

S^nt = np\na + 2n,„), K = /3/{2p{N/V)) 
(3 ^b/2 + 3Nf/4:-2, b = {ll/3)Nc - {2/3)Nf. (43) 

The free parameter k characterizes the diluteness of the ensemble. In the case under consideration with = 3 and 
Nf = 2 the value which reproduces the phenomenological diluteness of the instanton vacuum is k w 130 '18]. 

Within this assumption the color dependence in the activities only enters through the combination of the Tja^s 
which then can be integrated analytically, rendering the fermionic determinant real. For two flavor one obtains 



Za{zi,r) cx j dU Tia{T,p)T'\j^{T,-p) 



z^izi,r) cx j dU {Tia{T,p)T\^{T,-p))'^> 

^ (2iv. - i)(/+/r + f+f^f + {f+f^ - f^f+f 

AN^[N^ - 1) 

where ff^ = fi(zLp) are the functions defined in the previous section. The last free parameter of the model is the 
constant C that we fixed imposing that the absolute minimum of the thermodynamic potential Q{T = 0, p — O'.Ua, n^) 
appears at a total instanton density of N/V = Ha + 2n„i = 1.0 fm^* which is the phenomenological instanton density 
in the vacuum. 

The numerical results obtained for the Nc — 3, Nf = 2 cocktail model in the chiral limit, are presented in Fig. [5] 
The atomic density decreases as expected with growing temperature and chemical potential, while on the contrary 
the molecular component of the ensemble becomes more and more important. The behavior of the total density 
is represented in the upper panel: at small chemical potential when the temperature grows the instantons tend to 
diminish and eventually disappear at a temperature around T « 150 MeV. On the other hand, the presence of a 
non-vanishing density of quarks tend to preserve instantons. 

The mean field approximation allows establishing the connection between the density of the atomic component and 



the chiral quark condensate shown in ( 28 1 . From this latter relation one can connect restoration of chiral symmetry 
with the disappearance of atomic instantons. Indeed, our results, presented in Fig. [2] show a phase transition at 
p Ki 320 MeV for vanishing temperature and at T « 120 MeV at zero chemical potential. This numbers agrees with 
previous calculation in the instanton model [3 |T71 [13 [22] and confirm that the chiral transition appears when the 
total density of instantons is considerably different from zero in the full phase diagram. 



IV. CONCLUSIONS 



The instanton model has been extensively used in the study of the non-perturbative sector of QCD. The interest in 
it lies in the early discovery of the Dirac zero mode that the presence of instantons in the ground state of the theory 
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FIG. 2: Upper panel: Instanton density (red-yellow colors) and molecule density (blue colors) as a function of the temperature 
and chemical potential. Lower panel: total density (instantons + molecules). 

gives rise to. Despite all the development in the field, the general form of the zero mode profile in the general case of 
T ^ and /i ^ has been a pending issue. 

In this paper we fill this gap presenting the zero mode solution of the Dirac equation at finite temperature and 
chemical potential in the background of a singular gauge instanton. Starting from an Ansatz inspired by the known 
solutions in the limiting cases T 7^ 0, /i = and T = 0, /i 7^ we show that the proposed profile indeed satisfies the 
Dirac equation. 

The analytic expression for the zero-mode presented in this paper represents the starting point for the investigation 
of the role of instantons in the full phase diagram of QCD. In particular, it would be interesting to understand the 
contribution of such gauge configuration to the restoration of the chiral symmetry breaking around the critical point; 
on the other hand, also the U{1) axial anomaly and its possible restoration, can now be studied in the general setting 
of finite temperature and chemical potential. 

As a first step in this direction, we used the "cocktail" version of the model in order to investigate the instanton 
density in the phase diagram. In this framework, the analysis of the balance between the atomic and molecular com- 
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ponents of the ensemble provides some clues for unveiling the mechanism responsible for chiral symmetry restoration. 
In particular, we have shown that the restoration happens when the total density of instantons is still finite and, as 
expected, connected with the formation of instanton molecules in the vacuum. 

By studying the instanton model with more sophisticated techniques than the simple "cocktail" approach used here 
(e.g. the Interacting Instanton Liquid model developed in [7 ), we can hope to extract more information on the chiral 
dynamics induced by instantons in the QCD phase diagram. 
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Appendix A: Properies of 77 symbol 

The 't-Hooft 77-symbols entering in the analytic expression of the instanton profile, are defined through [5] 

The ?7-symbols are (anti) self-dual in the vector indices 

Contraction involving the 77-symbols can be done using the following relations: 

VafJ't^VapX — ^j.Lp^uX ~ ^^X^vp ~t~ ^^upX-, 
Vaj-^J^Vhj-ip — ^ab^up ~t~ ^ahc^cy pi 

^ abcVbfiiyVcpX ^ ^ fip^aiyX ^ pX^avp ~^ ^I'X^afj.p ^up'^apX-, 
^ XjAUcrVapa ^ pXVafii/ ^ pu^aXp ^ ppVah'Xi 

with the same relations holding for r/^^ix, except for 

VapvVapX ^^p^uX ^ ^pX^vp ~ ^jAUpX- 



Appendix B: Gamma and sigma matrices 

Let us define 

7f=7f, 7,f.--7^^ m = l,2,3. (Bl) 
where the explicit form of the gamma matrices is given by 



7a. 



-lal 
icr^ 



with = ((?, —i) and cr's the usual Pauli matrices with = 6°"^ + i(°'^'^a'^. The resulting 75 will be 



75 = 71727374 = 



1 
-1 



Directly from the definition we have 
so that 

Notice finally that the 't Ifooft 77 symbols can be defined using the cr^ by 



1 
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